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Abstract. - Almost identical thermal properties of type-I clathrate compounds to those of glasses 
follow naturally from the consideration that off-centered guest ions possess electric dipole moments. 
Local fields from neighbor dipoles create many potential minima in the configuration space. A 
theoretical analysis based on two-level tunneling states demonstrates that interacting dipoles are 
a key to quantitatively explain the glass-like behaviors of low-temperature thermal properties of 
type-I clathrate compounds with off-centered guest ions. 
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Introduction. — Clathrate compounds, attracting 
much attention in connection with thermoelectric mate- 
rials [1], show unique physical properties at low tempera- 
tures. A series of recent experiments on type-I clathrate 
compounds, /3-Ba 8 Ga 16 X 30 (X=Sn,Ge) or /J-XgGaigGeao 
(X=Eu, Sr, Ba) — hereafter abbreviated by /3-BGX, (3- 
XGG, or /3-BGS, etc. — have revealed that the character- 
istics of thermal [2-6] and spectroscopic [7,8] properties at 
low temperatures strongly depend on the states of guest 
ions in cages. These studies confirmed earlier measure- 
ments on thermal properties for various types of clathrate 
compounds [9-15]. 

The most intriguing findings [2-15] are that the specific 
heats and the thermal conductivities of type-I clathrate 
compounds with off- centered guest ions show almost iden- 
tical thermal properties to those of glasses in the regime 
below a few K [16], in spite of the remarkable dif- 
ference in microscopic structures among these materi- 
als. The specific heat of /3-BGS is scaled by the rela- 
tion C ~ aT + /3T 3 with a ~ 30 mJmol^K- 2 and 
(3 ~ 50 mJmol _1 K -4 [2-6]. While topological disorder 
is a key element for glasses [17], the disorder in type-I 
clathrate compounds is much less relevant, as diffraction 
studies have provided the evidence that the microscopic 
structures of type-I clathrate compounds are rather well- 
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defined [2-7,11,13]. 

A feature found in experiments [2-15] is that guest ions 
take either an on-center or an off-center position depend- 
ing on the size of the cages or, equivalently, on the ionic 
radii of the guest ions. It is notable that the thermal 
properties almost identical to those of glasses are only 
observed in clathrate compounds with off-centered guest 
ions. It remains a mystery why the observations for these 
clathrate compounds without topological disorder show 
thermal properties almost identical to those of glasses [16]. 

The experimental data so far have been analyzed on 
the basis of the model that the isolated off-centered guest 
ion governs the thermal properties at low temperatures. 
This simple view often misleads us to the idea that a 
noninteracting guest-ion picture could explain the glass- 
like behaviors. Such a noninteracting guest-ion picture, 
however, conflicts with the observation, yielding a specific 
heat two orders larger than the observed value at 1 K, as 
demonstrated below. This Letter provides a clear inter- 
pretation of the glass-like specific heat behaviors of the 
type-I clathrate compounds with off-centered guest ions, 
/3-BGS [2-6]. Our treatment is general and directly ap- 
plicable to other types of clathrate compounds with off- 
centered guest ions. 

Off-centered guest ions and tunneling states. 

Experimental studies have shown that the positional sym- 
metry of the guest ions in clathrates is broken by increas- 
ing the size of the cages, and it has been experimentally 
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confirmed that the guest ions take four off-center positions 
at ro from the center of a cage [2-7, 11, 13]. A typical dis- 
tortion ?'o = 0.43 A in /3-BGS is obtained from diffraction 
experiments [2,4]. This value is about 7.2% of the dis- 
tance between the neighboring 14-hcdrons (a/2 = 5.84 
A, where a is the lattice constant). It is straightforward 
to see that such an off-centered guest ion induces a large 
electric dipole moment due to the difference of the charges 
between the Ba 2+ guest ion and the Ga~ ion constituting 
cages. The strength of the electric dipole moment be- 
comes p w 4.1 Debyes in /3-BGS from ro = 0.43 A. Here, 
we emphasize that the cages with off-centered guest ions 
intrinsically possess dipoles and that it is crucial to prop- 
erly deal with their characters in modeling the systems. 

It has been also confirmed [13] that guest ions experi- 
ence a hindering potential Vh{9) with a four- fold inver- 
sion symmetry along the azimuthal direction. The bar- 
rier height of the hindering potential Vh between nearby 
potential wells is estimated from the difference between 
the largest and the smallest strengths of van dcr Waals 
potentials to a guest ion from the framework atoms of a 
cage, which should be of the order of 10 K. Actually, first- 
principles calculations have shown Vh to be ~ 20 K for 
Sr 2+ guest ions in /3-SGG [18]. By expressing the posi- 
tion of a guest ion in terms of two-dimensional spherical 
coordinates perpendicular to the four-fold inversion axis 
centered in a cage, the motion of an off-centered guest ion 
is mapped onto a rotational motion of dipoles having a mo- 
ment of inertia I = 25.4 uA 2 in the case of j3— BGS, where 
the kinetic energy is scaled by Ek = h 2 /(2I) = Ek = 0.96 
K. Thus, we have the situation Ek < Vh- 

Here, we mention that /3-BGG are slightly off-centered 
by 0.15 A [19, 20] with much smaller dipole moments 
compared with the case of /3-BGS. In this connection, it 
is noteworthy that p-type /3-BGG shows crystalline-like 
thermal conductivity, while n-type /3-BGS behaves glass- 
like [2] . This is one of the evidences that the long-range 
dipolc-intcraction is crucial for interpreting the thermal 
properties of type-I clathrates since the Coulomb inter- 
action between the cages and guest cations Ba 2+ in (3- 
BGG is shielded in n-type electron-rich BGG, while the 
dipole interaction in p-type becomes relevant even for 
small dipole moments. The network configuration con- 
sisting of the off-centered guest ions is schematically illus- 
trated in Fig. 1. 

The separation of neighboring wells is small, for exam- 
ple, 7rr /2 = 0.67 A in /3-BGS. This allows the guest 
ion to tunnel to a nearby potential minimum at lower 
temperatures than the barrier height energy Vh, where 
the off-centered guest ion executes zero-point motion at 
one of the four wells with an energy huiQ. We can take 
into account only two wells with a small energy difference 
e = Ei — E2 (the asymmetric energy of neighboring double 
wells), which mainly attribute to the disordered configu- 
ration of the Ga~ and Sn framework atoms of the cages. 

In this case, the Hamiltonian is reduced to a 2x2 ma- 
trix with the off-diagonal element A (the tunneling term 
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Fig. 1: (Color online) Schematic illustration of the configura- 
tion of the guest ions in the over-sized cages in /3-BGS. The 
four-fold inversion axes are directed along x, y, z. (a) The de- 
viation from the on-center position (filled circles) induces the 
electric dipoles (arrows). Here the two nearest dipoles are de- 
picted. The electric dipoles rotate in the plane perpendicular 
to the axis linking the nearest dipoles. (b) The configuration 
of the on-center positions of the guest ions. The filled circles 
represent the positions of off-centered guest ions, around which 
the electric dipoles are induced. The sites A, B, and C in (b) 
are seated on the chains parallel to x, y, and z, respectively: 
A = (a/4, 0, a/4), B = (0, a/4, 3a/4), C = (a/2, a/2, a/2). The 
distance between the next-nearest neighbors (dashed lines) is 
-\/3/8 a. Note that these three dipoles constitute an equilateral 
triangle and easily generate a frustrated situation, (c) The 3D 
configuration of the dipole chains is illustrated. 



due to the overlap of the wavefunctions), and the ener- 
gies are given by E = ±(l/2)Ve 2 + A 2 . The specific heat 
incorporating all of the off-centered guest ions is given by 



C 



2k%T 



x 2 sech 2 {x)P{2xk B T)dx, 



(1) 



where x is the dimensionless variable defined by 
-E/(2fcsT), and P(E) is the number density of states per 
volume; namely, P(E)dE is the number of tunneling states 
per volume with energy between E and E + dE as de- 
scribed in [21,22]. 

Noninteracting-dipole picture. The picture that 
tunneling between nearby sites in Vh(0) is relevant to the 
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glass-like thermal properties often leads to the erroneous 
idea that the noninteracting-dipole picture would work 
to explain the glass-like behaviors. Such a model, how- 
ever, does not reproduce the observed temperature depen- 
dence and the magnitudes of the specific heats for typc-I 
clathratc compounds with off-centered guest ions as shown 
below. In fact, as will be discussed in the next section, the 
tunneling between the nearby potential minima in the con- 
figuration space generated by interacting dipoles is crucial 
for interpreting glass-like thermal properties below about 
1 K . 

To disconfirm the noninteracting-dipole picture, we cal- 
culate the specific heat within this picture. When the 
dipole interaction between the guest-ions is irrelevant, the 
problem reduces to the state of only one cage. In this 
case, the zero-point energy of the guest ion trapped in 
one of the four local-potential minima of the hindering 
potential Vh(0) is estimated to be Tiuiq r* 4.6 K from 
the uncertainty principle, taking an actual size of box 
~ (0.3 A) 3 for /3-BGS [18]. The tunneling energy is given 
by A = huj exp(—y/2IVh.S6/h) with 66 the angle between 
two nearby minima, and 69 takes a value ~ n/2. If the 
asymmetry energy of neighboring wells e is negligibly small 
compared with A, we can estimate for the noninteracting- 
dipole picture the most probable lower-bound energy as 
A m i„ = 0.03 K by using values for /3-BGS: huiQ = 4.6 K 
and 69 = ir/2. Note that, in the case of the interacting- 
dipole picture, A m i n becomes much smaller (see below). 

The lower and the upper bounds of the integral in 
Eq. (1) should be A mm /(2k B T) and E max /(2k B T), re- 
spectively, depending on the temperature T. Since the 
function x 2 sech 2 (x) in the integrand has a maximum at 
around unity, the contribution to the integral of the prod- 
uct of x 2 sech 2 (x) and P(2xk B T) becomes small at tem- 
peratures below T < A m i n . In addition, the noninteract- 
ing-dipole picture claims that every off-centered guest ion 
contributes to the specific heat. By taking the width E of 
the P(E) from A m i n to E max w 2Vh and by incorporating 
the number of guest ions in a unit cell 6, Eq. (1) yields 
C ~ 4100 mJ mol -1 K _1 at 1 K, which is two orders larger 
than the observed value for /3-BGS C ~ 30 mJmol -1 K _1 
at 1 K [2,4,5]. Thus, the noninteracting-dipole picture 
severely conflicts with the observations. 

In this connection, it is noteworthy that the Einstein 
model has been used to analyze the Boson peak-like be- 
havior of specific heats or inelastic neutron scattering data 
observed at the around 10 K energy region. As is well 
known, the Einstein model is a good phenomcnological 
approximation for describing the excess phonon density of 
states, but is not suitable to grasp the essential physical 
picture behind the model such as the symmetry or many- 
body interaction involved in the systems. 

Interacting-dipole picture. — Electric dipole mo- 
ments of guest ions provide long-range interactions with 
dipoles in other cages. To make our argument clearer, let 
us consider at first two electric dipoles pi and pj sepa- 



rated by a distance \Rij\- The dipolar interaction gives 
the following form under the condition rg <C Rij: 



VijiPhPj) = 



Pi ■ pj - 3(pi ■ //, /.■■://, ■ Rjj) 



(2) 



where e r is the dielectric constant of the clathrate con- 
sidered here, and Rij = Rij/\Rij\. In the case of the 
/3-BGS clathrate, off-centered guest ions are in 14-hedron 
cages, and its nearest-neighbor off-centered guest ions are 
located in the next two 14-hedron cages, which share the 
same four-fold inversion axis. Figure 1 gives the configu- 
ration of the guest ions in over-sized cages. The four-fold 
inversion axes are directed along the x, y, z axes due to the 
cubic symmetry of /3-BGS. 

The potential function for the two coupled dipoles p\ 
and p 2 becomes V 12 = V h {6{) + V h {9 2 ) + Wi 2 (9 1 ,9 2 ) with 
W12 = p 2 cos(#i — # 2 )/(47T£vi?f 2 ), where two global min- 
ima (maxima) in (9i,9 2 ) configuration space appear at 
\9i — 9 2 \ = n (2tt) since the dominant term for a nearest- 
neighbor pair is the first term in Eq. (2). This configura- 
tion acts as a new hindering potential in addition to the 
four- fold inversion symmetric potential Vh- This argument 
can be straightforwardly extended to the case of multiple 
pairs providing many potential minima in the configura- 
tion space V = (6*i, 9 2 , 6*3, • • •) [23], where the potential 
function is V 123 ... = EW) + E 

The energy scale of the dipolar interactions between 
nearest neighbors is given by its maximum value: 



■h 



P 



4TTe r Rf 



(3) 



where i?i is the distance between the nearest neighbors. 
We also represent the energy scale of the next-nearest 
neighbor coupling as J 2 ~ Ji/2. Here, we defined J 2 as 
the maximum value of the coupling energy averaged over 
the eight next-nearest neighbors. The contribution to J 2 
from the first term in Eq. (2) is smaller than that of the 
second one, and mostly works as a counterforce. Note 
that the motion of the dipoles is restricted in the plane 
perpendicular to the nearest-neighbor chain. Therefore, 
wc first estimate J 2 given by J 2 = 3Ji{Ri/R 2 ) 3 (()max 
without the first term of Eq. (2), and then take it into 
account, where R 2 (= \/3/8 a) is the distance between 
the next-nearest neighbors, and (Q m ax the maximum av- 
erage of C = cos(#i) cos(#2) over the angles. Setting 9 2 as 
a function of Q\ which yields the largest ( for each 9\ , we 
average £ over 6\. Since 9\ takes a value between 9 m i n 
and (ir — 6 m i n ), where 9 m i n = tan _1 (l/v / 5) in /3-BGS, it 
follows (Qmax — 0-5 and J 2 ~ O.8J1. The average of the 
first term over eight neighbors becomes ~ — 0.3Ji in the 
case that four next-nearest dipoles are anti-parallel and 
the others are perpendicular to the central dipole. Since 
the contribution to J 2 from the first term is not far from 
this, we can take J 2 ~ 0.5Ji. 

The actual distance between nearest-neighbor guest ions 
in /3-BGS is i£i = a/2 = 5.84 A, then the characteristic 
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energy scale for nearest-neighbor coupling is estimated as 
Ji ~ 6e r /eo K. Taking into account the dielectric con- 
stant for semiconductors in the range 5 < e r /eo < 20, 
it turns out that the guest ions experience strong interac- 
tions. They are no longer regarded as independent dipoles. 
Combining four-fold inversion symmetry of dipoles with 
the frustrated situation due to the equilateral triangle (Sec 
Fig. 1), many local minima are created in a hierarchical 
potential map in a configuration space V, where tunnel- 
ing should involve simultaneous local structural rearrange- 
ments of an appropriate number of guest ions. 

An EXAFS study [19] for /3-BGG has shown that the 
off-center sites actually have a three dimensional (3D) 
components. Though our analysis is based on the 2D 
picture, our main conclusions are not altered even taking 
account of 3D character since the key aspect is to gener- 
ate multi-valley potentials in the configuration space by 
long-range dipolc-intcraction. 

P(E) should be symmetric about E = since E\ > E% 
or Ei < Ei is equally likely, implying P(E) is an even 
function of E. It is reasonable, at the first glance, to 
take P{E) as a normal distribution at the center of the 
distribution. We employ the form of P(E) as 

P(£) = J^exp[-(£/E) 2 ], (4) 

where V = Nacl 3 (Na- Avogadro number), and n is the 
number of relevant tunneling states per volume, which we 
set to n = 677/a 3 by incorporating the number of off- 
centered guest ions, 6, in a unit cell a 3 . The numerical 
factor 77 expresses how many of the off-centered guest ions 
are relevant to the specific heat in effect: Note that 1/?/ 
is the averaged number of off-centered guest ions involved 
in a rearrangement of the configuration. 

Explicit form of T-linear specific heat. By 

taking the upper bound (~ 2E) from Eq. (4) as the 
maximum coupling strength from neighbor dipoles 2E ~ 
z\J\ + Z2J2 — 6J1 with configuration numbers Z\ = 2 and 
z-i = 8 (J3 <C J\, J2), and noting n/J\ = 2>ite r ri/p 2 from 
R\ = a/2, Eq. (1) combined with Eq. (4) yields the simple 
form 

for ksT > Amin- Equation (5) predicts larger specific 
heats for smaller dipole moments p. Suekuni et al. [5] have 
observed low-temperature specific heats for n- and p-type 
/3-BGS with different length dipoles of ro = 0.434 and 
0.439 A, respectively. The observed specific heats below 1 
K of the n-type /3-BGS are a few percent larger than those 
of the p-type /3-BGS. This accords with the prediction of 
Eq. (5). ' 

The comparison of Eq. (5) with the observed magnitude 
of a specific heat C ~ 30T mJ mol" 1 K" 1 for /3-BGS [2-5] 
below 1 K [2-6] provides important information for tun- 
neling states through the value of rj. Taking e r — 10eo, 



r] is estimated to be ~ 0.06 from C ~ 50(e r /eo)r]T 
mJrriol -1 K _1 . This indicates that not all but only 6% 
of the off-centered guest ions contribute to the specific 
heats on average. This implies that the averaged number 
of dipoles simultaneously rearranged is l/rj ~ 20, which 
indicates that A m i„ is negligible in the temperature range 
we consider here. 

The glass transition temperature T g should be close to 
T g ks 2T me i t /3 f=a 500 K from the Tammann's rule, where 
Tmeit ~ 4E(~ 700 K) estimated from the effective width 
of the distribution function Eq. (4). This temperature is 
much lower than the one where the dipole moments dis- 
appear in /3-EGG (> 1000 K) [7]. We suggest an exper- 
imental study to confirm the existence of glass transition 
in this temperature range. 

Summary. — The results of our investigation high- 
light the role of the off-centered guest ions in /3-BGS, pro- 
viding a clear picture of why this compound shows glass- 
like behaviors. The off-centered guest ions constitute tun- 
neling states in a configuration space of interacting dipoles 
with hierarchial structure, where tunneling occurs as a re- 
arrangement of the dipole configuration. The rearrange- 
ment is induced in units of 5-20 dipoles. Although these 
results are obtained on /3-BGS, the investigations are di- 
rectly applicable to other types of clathrate compounds 
with off-centered guest ions. 

We have not discussed about the glass-like behaviors of 
thermal conductivities observed in type-I clathrates on the 
basis of our picture. This subject will be treated in detail 
in our forthcoming paper. 
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